Griffiths phase in diluted magnetic semiconductors 
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We study the effects of disorder in the vicinity of the ferromagnetic transition in a diluted magnetic semicon- 
ductor in the strongly localized regime. We derive an effective polaron Hamiltonian, which leads to the Griffiths 
phase above the ferromagnetic transition point. The Griffiths-McCoy effects yield non-perturbative contribu- 
tions to the dynamic susceptibility. We explicitly derive the long-time susceptibility, which has a pseudo-scaling 
form, with the dynamic critical exponent being expressed through the percolation indices. 
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Introduction. There has been increasing interest recently 
in studying diluted magnetic semiconductors both experimen- 
tally and theoretically 0|. This interest is motivated not 
only by the possibility of potential technological applications 
of the materials but also by a very rich physics, which follows 
from the unique combination of magnetic and semiconduc- 
tor properties co-existing in diluted magnetic semiconductors. 
There is a large number of different materials currently be- 
ing studied such as Gai_ x Mn x As, Ini_ x Mn x As, Gai_ A Mn x P, 
Gai_ x Mn c N, Ge\- x Mn x , etc. One of the common important 
features of these materials is the invariable presence of strong 
disorder, which plays an essential role in both magnetic and 
transport properties of the systems with the most prominent 
effect being the localization of carriers. Disorder arises here 
from the random locations of magnetic impurities in the host 
lattice of the semiconductor and also from other impurities 
and defects in the system. It is important to emphasize that 
disorder effects are completely neglected in the continuum vir- 
tual crystal approximation often used in the theoretical studies 
of these systems. 

It is now widely accepted that ferromagnetism in diluted 
magnetic semiconductors is due to an indirect interaction be- 
tween magnetic impurities, which is mediated by holes. The 
effective magnetic coupling strongly depends upon the ratio 
nV 3 Li oc of the localization radius and the mean distance be- 
tween the magnetic ions. In the strongly localized regime, 
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i.e. when n. ' L\ oc <sc 1, the system may be viewed as a com- 
bination of two static magnetic subsystems (magnetic impu- 
rities and holes) with the corresponding Hamiltonian having 
the following form: 



'H = Y J J d r « - «■«!) S««Tj +Yj ^af (|r fl - r h \) S a S h , 



(1) 



where S fl are impurity spins and er, are hole spins, J{r) is the 
manganese-hole interaction, which is given by 



J(r) = /o exp 



-4oc 



(2) 



and Jaf(>") is the direct antiferromagnetic coupling originating 
from the exchange interaction of the deep-lying electron states 
of magnetic impurities, which decays exponentially with dis- 
tance much more rapidly than the manganese-hole coupling 



(0. The Hamiltonian Q leads to a polaron percolation pic- 
ture |2B, which qualitatively explains the mechanism of the 
ferromagnetic transition in the system: In diluted magnetic 
semiconductors, the concentration of holes is much smaller 
than the one of manganese atoms If T <K \J$\, the magnetic 
impurities around each localized hole get polarized forming 
a bound magnetic polaron. As we lower the temperature, the 
polaron size R p (T) increases and polarons overlap forming fi- 
nite clusters. At some threshold temperature T c , an infinite 
cluster appears signaling the ferromagnetic transition. The di- 
rect antiferromagnetic coupling cuts off the polaron growth at 
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the distance R* = Li oc In |7o/-/af("- )|, when the manganese 
impurities prefer the antiferromagnetic alignment to the cou- 
pling to a distant hole. If this distance exceeds the polaron size 
at the transition point R* » R p (T c ), the antiferromagnetic in- 
teraction does not affect the ferromagnetic transition physics 
and can be neglected at T ~ T c (we accept this assumption in 
the present paper). 

It is possible to formalize the polaron picture developed in 
the paper 1 3] by deriving the effective Hamiltonian for inter- 
acting polarons. This can be achieved (see below) by integrat- 
ing out the manganese spins. The resulting Hamiltonian (jSJl 
is the ferromagnetic Heisenberg-type Hamiltonian with small 
logarithmic "non-Heisenberg" corrections, which can be ne- 
glected. The corresponding interaction constant depends upon 
the inter-polaron distance, which is the distance between the 
localization centers of the corresponding holes. Since the 
spatial distribution of the holes is the random Poisson dis- 
tribution, the ferromagnetic polaron coupling becomes a ran- 
dom variable with a well-defined distribution function. This 
maps the initial problem onto a strongly disordered ferromag- 
netic Heisenberg model with randomly placed sites. This sys- 
tem naturally leads to the well-known Griffiths-McCoy effects 
0, El in the disordered phase near the phase transition. The 
Griffiths-McCoy singularities are central of the present Letter. 

The original paper of Griffiths |4] addressed the nature of 
the transition in a random Ising ferromagnet, although the 
main conclusions are relevant to much more general systems. 
The central idea of the Griffiths theory is that above the fer- 
romagnetic transition point in a disordered system, there is 
always a finite probability of finding an arbitrary large ferro- 
magnetic cluster. These rare-fluctuation clusters give singular 
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contributions to the magnetization, which can be proven to be 
a non-analytic function of the external magnetic field. The 
corresponding phase is called the Griffiths phase. The state- 
ment of the non-analyticity of the magnetization is a purely 
mathematical statement. The corresponding corrections, al- 
though singular, are extremely (exponentially) small and it 
is not clear whether such effects in thermodynamic quanti- 
ties are observable in real experiments. However, the Grif- 
fiths effects are much more pronounced in dynamic quantities 
because the disorder in the model is quenched (does not de- 
pend on time) and thus the rare-fluctuation induced clusters 
are infinitely extended in the time direction. Thus, in the Grif- 
fiths phase region, the dynamic response is governed by the 
Griffiths-McCoy effects as was shown in Ref. [6]. Currently 
there is a great deal of interest in Griffiths physics related to 
quantum phase transitions 01 (a quantum cZ-dimensional 
system is equivalent to a d + 1 -dimensional classical systems 
at T — 0; thus, a quantum phase transition in a disordered sys- 
tem is governed by Griffiths-McCoy effects due to disorder 
being perfectly correlated in time). 

Effective polaron Hamiltonian. Let us consider two local- 
ized holes embedded into the system of impurity spins. Let the 

— 1/3 

inter-hole distance be large: r » n. . The partition function 
for the system reads: 



Z = Tr trij0 . 2 Tr {Sa} | _ [exp|--!§ a [7(|r a -ri|)d-i (3) 



+J(\r a -r 2 |)<x 2 ] 



It is important to realize that we are building an effective de- 
scription in terms of bound magnetic polarons, which are very 
"heavy" classical objects, and the impurity spins are generally 
large (5/2 for manganese impurities). Thus, we can ignore 
quantum effects and replace the traces in Eq. by the corre- 
sponding classical integrals. We rewrite Eq. as 



Z = 



9 



An An 



exp<- 



£ e ff (cos 6, T) 



(4) 



where integration over £1, implies averaging over the orienta- 
tion of the ith polaron's spin, 9 is the angle between the spins 
of the two polarons, and the effective interaction is defined as 
[see Eq. J3J]: 



fieff ( cos S, T) 



T y ln ( sinh[7(r fl )/r] 



(5) 



with 



J(r) = Jo V2exp 



r\ + r 2 1 



2L 



jcosh 



loc 



n - r 2 



+ cos 9, (6) 



where ri >2 are the distances between a magnetic impurity and 
polaron centers (localized holes). Taking into account the fact 
that the density of magnetic impurities is high, we can replace 



the sum in Eq. (|5J by the integral: 



<9£ e ff (cos 9, T) 
d cos 9 



= -m j d i 



r tanh 



J(r a ) 



dJ(r) 

dcosO 



(7) 



In the limit L\ oc <K r, we have 



£ eff (cos0,r)= - ^g4 c rni)exp(-^- 



xln 



f(6)je^ 



cos 9, 



(8) 



where f(9) ~ 1 is a function of 9, which is always non-zero 
and can be neglected within the logarithmic accuracy. We 
conclude that the effective polaron Hamiltonian is the ferro- 
magnetic exchange Hamiltonian with randomly placed cites 
and the couplings between them exponentially decaying with 
distance. 

A similar random exchange Hamiltonian was considered by 
Korenblit et al. in Ref. |9], where it was found that the fer- 
romagnetic transition in the system can be conveniently de- 
scribed with the help of the percolation theory. Within the 
percolation approach, two interacting polarons are considered 
either locked in the same direction (if the effective interaction 
is larger than temperature: > T) or completely uncorre- 
cted in the opposite limit, < T (which essentially means 
that the weak couplings can be set to zero). This means that 
we replace the initial system of randomly interacting polarons 
by a system in which the weaker couplings are cut. When 
decreasing temperature, the number of strongly coupled po- 
larons increases and at some temperature, an infinite cluster 
appears. This temperature is identified with the ferromagnetic 
transition temperature of the system. We would like to em- 
phasize that the percolation theory is able to establish T c only 
up to a numeric factor of the order of unity. In the immediate 
vicinity of the ferromagnetic transition point, the percolation 
theory may fail since it is not able to account for critical fluc- 
tuations, which become increasingly important at T — > T c . 
However, there is a small parameter in our theory, namely 
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n h L] oc <k 1. The parameter p of the percolation theory is 
related to temperature by 



^)^< /3 ilodn(y) , 



(9) 



with p c ss 0.86. Thus there exists a domain 1 » \p— p c \/p c » 
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" h L\oc, which is outside the thermodynamic critical region 
\T - T c \/T c <s 1, but still close to the percolation transition, 
where we can use the asymptotic results of the percolation 
theory. 

Griffiths theorem. Let us prove the existence of the Grif- 
fiths phase in the framework of the percolation theory. In the 
percolation picture, we have a number of finite clusters above 
the transition point. The concentration of clusters containing 
N polarons for N » % D is iMIuil 



P(N, T) cc N- T ^ | exp 



N 



T-9 



(10) 
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where N is the number of polarons in a cluster, ^ oc 
(\p(T) - p c \ lp c ) v is the dimensionless correlation length. 
Percolation exponents are as follows 

r * 2.18, 0.88, (11a) 
T > T c : £ = 1, 6 = 3/2, D = 2, (lib) 
T < T c : £ = 2/3, 6 = -1/9, D = 3 (11c) 

for the three-dimensional percolation. Constant A in Eq. (II Oi 
is a number of the order of unity. 

Each polaron has a spin S p — S n\v p , where S is the spin of a 
single magnetic impurity and v p is the polaron volume which 
is the volume of the sphere of radius R p L\ oc ln(Jo/T). The 
magnetization of the system in an external magnetic field H 
can be written as follows: 



M(H) = S P Y J NP(N) 



coth 



NS„H 



NS P H 



Let us expand this function in Taylor series 



(12) 



(13) 



Following Ref. (Till , we obtain the radius of convergence R H 
of the series: 





lim 






Up, 


N— >+oo 


on 
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Above the percolation threshold, ( = 1 and the series Jl 3i has 
zero radius of convergence, which implies the non-analyticity 
of the magnetization as a function of the external magnetic 
field due to rare disorder fluctuations, i.e. the essential Grif- 
fiths singularity. The corresponding singularity in the magne- 
tization and susceptibility is however very weak: 



M(H) oc exp 



-const 



(!) 



(14) 



which is typical for the Griffiths effects on thermodynamic 
quantities in classical systems. 

Dynamic susceptibility. The Griffiths physics manifests it- 
self much stronger in dynamic quantities. The quantity of in- 
terest is the dynamic susceptibility, which is connected with 
the spin-spin autocorrelation function 



C(?) = ^:Z <SflWSfl(0)> ' 



(15) 



In the Griffiths phasejjhe disorder-averaged correlator can be 
written as follows I 



C(r) = J> (A ° Cw(f) ' 



(16) 



where Cj</(t) is a function describing the relaxation of the mag- 
netic moment of a cluster of size N. 



The relaxation of the magnetic moment in a cluster is gov- 
erned by the Gilbert equation: 



dm 
~~dt 



ym x 



dm 



dm 



(17) 



where y and r\ are constants corresponding to the gyromag- 
netic ratio and Gilbert damping and m = M/Vjv, is the mag- 
netization density with Vn being the volume of the cluster 
V N = Nv p = A^(4^/3)Lf oc ln 3 (y /r). Function T{m) is the 
free energy of the cluster per volume as a function of the di- 
rection of the magnetic moment, and h(f) is the fluctuating 
magnetic field. Let us note that the case of h — and !F(m) — 
-Hm yields the well-known Landau-Lifshitz equations fol 
with the dimensionless damping coefficient a = rjym. We 
assume that anisotropy of the cluster is weak enough so that 
parameter [Vn {Tmax — ^min) /T] is very small, which corre- 
sponds to the superparamagnetic limit. 

The random magnetic field h(f) is due to thermal fluctua- 
tions. They are supposed to be Gaussian, with the following 
white-noise correlation functions: 

<h(f)> = and (KWpiO)) = p6 ap 5(t), (18) 

where Greek indices label Cartesian coordinates. It is possible 
(see Ref. f3l) t0 connect the constant fi with the properties 
of the cluster and temperature. To do so, one imposes the 
natural condition that the stationary solution of the stochastic 
Gilbert equation fl!7i subject to the thermal noise J 1 81 reduces 
to the Boltzmann equilibrium distribution for the magnetiza- 
tion density: m(9, </>) oc exp[-y/ v i^ r (6', </>)/r]. This leads to the 
following expression for the constant in the correlator Jl 81 : 
H = 2i]T/V N . 

It is further possible to derive general Fokker-Planck equa- 
tions for the distribution function of the magnetic moment 
P[m; t] in the case of an arbitrary anisotropy ^"(m) (this was 
done in the pioneering paper of Brown I13I0 . We however 
confine ourself to studying the isotropic case only when the 
Fokker-Planck equation reduces to the simple diffusion equa- 
tion on a sphere (also described phenomenologically by Bray 
in Ref. 1 6]): 



P[m, t] = 0, 



(19) 



where Ag t <f, is the angular part of the Laplacian and the effec- 
tive diffusion coefficient is equal to (see e.g. Ref. lfl3ln 

mT) --U{rh)^^ (20) 

T 1 loc 

where a is the dimensionless damping coefficient introduced 
in the previous paragraph, y is the gyromagnetic ratio, and S 
is the magnetic impurity spin. 

Let us suppose that the initial magnetization m(t = 0) of the 
cluster is in the z-direction, which corresponds to the initial 
condition P[9, (/>; t = 0] = 6(6) for the distribution function, 
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defining the Green function for the diffusion equation il9\ . 
The corresponding solution reads: 

P[0\t] = ^]Pi (cos 8) exp [-£>/(/ + l)f] sin 0, (21) 

where f/( ) are Legendre polynomials. The z-component of 
the magnetization averaged over the distribution function (12 1 i 
is exactly the function C^(f) introduced in Eq. dl6> . which 
describes the relaxation of the magnetic moment in a cluster 
of size A^: 



C N (t) = Mexp [-2D(N, T)t] . 



(22) 



We now proceed to calculate the spin autocorrelation func- 
tion by replacing the sum in Eq. i ]16> by the integral 



C(t) 



r n [ - 6 

c aW— — - 



exp 



2D(N, T)t 



(23) 



which can be evaluated by the saddle point approximation: 



C(f) oc 



/ P(r)-Pc 



f fe exp 



, (24) 



where, for the sake of brevity, we have introduced the follow- 
ing indices: 



c 



b = 



1+f 
2-6- £/2 



1+f 



1 +f \2 

and the function in the exponent reads: 



g(r)=AT 



[2ND(N,T)]~> 



(25a) 
(25b) 
(25c) 
(25d) 

(26) 



where the diffusion coefficient is defined by Eq. (I20> and A is 
the number of the order of unity introduced in Eq. JlOl i. The 
explicit values for the parameters above (below) the transition 
are: A = 1/2 (2/5), c * 0.88 (1.06), b * (14/15), and 
a * 0.32 (3.29). 

Conclusion. We have shown that a diluted magnetic semi- 
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conductor in the strongly localized regime {i.e., if n. L\ oc «: 



1) must be in the Griffiths phase in the vicinity the ferromag- 
netic transition point. The magnetization is a non-analytic 
function of the external magnetic field with the correspond- 
ing singularity being exponentially small. The dynamic re- 
sponse in the vicinity of the transition point is determined by 
rare-fluctuation Griffiths droplets. These droplets correspond 
to the percolating clusters in the polaron percolation picture. 
The dynamic susceptibility has the form J24I of the stretched 
exponential In C(t) oc -t A , with A being connected to the per- 
colation indices and equal to 1/2 in the Griffiths phase and 
2/5 in the ferromagnetic phase. An important consequence 
of our work relates to the role of disorder in ferromagnetic 
semiconductors. Since disorder is fundamental to the occur- 
rence of Griffiths phase, an observation of Griffiths physics in 
diluted magnetic semiconductors will have significance with 
respect to the role of quenched disorder in these systems. In 
fact, we believe that insulating diluted magnetic semiconduc- 
tors {e.g., GaMnAs or InMnAs on the insulating side of the 
metal-insulator transition) may be an ideal system to look for 
signatures of evasive Griffiths physics. 
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